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Preface 


The  work  described  in  this  thesis  accomplishes  two 
tasks.  First,  two  electrocardiograms  (ECG)  of  a heart 
progressing  from  near  normal  to  P wave  disassociation  were 
recorded.  This  furnished  the  data  to  test  an  algorithm 
designed  to  identify  P waves.  The  algorithm  relied  on  the 
P wave  being  associated  with  the  R wave  only  during  the  first 
two  heart  beats.  After  this  a routine  of  "stepping"  the 
distance  between  P waves  (P-P  interval)  down  the  ECG  was 
used.  After  each  "step"  a search  was  made  for  the  peak  of 
the  P wave. 
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Abstract 

Heart  failure  was  induced  in  one  dog  by  a blockage  of 
the  anterior  descending  branch  of  the  coronary  ostium. 

Heart  failure  was  induced  in  another  dog  by  overdose  of 
digitalis.  The  electrocardiogram  {ECG)  from  each  dog  was 
recorded  on  an  analog  magnetic  tape  using  an  AMPEX  Model 
FR1300  recorder.  An  algorithm  to  locate  and  identify  the 
P wave  using  the  periodicity  characteristic  of  the  P wave 
was  tested  against  this  data.  The  P waves  in  the  first  two 
heart  beats  are  identified  by  their  proximity  to  the  R wave. 
The  distance  between  these  wave  shapes  CP-P  interval)  is 
"stepped"  down  the  ECG.  At  each  "step"  a search  is  initiated 
for  the  peak  of  the  P wave.  If  a peak  is  not  found  within  a 
set  number  of  data  points,  the  last  P-P  interval  is  retained. 
The  accuracy  of  the  algorithm  is  excellent  until  an  extreme 
abnormal  wave  shape  interferes  with  the  ECG  around  the  area 
of  the  P wave.  When  this  happens,  the  algorithm  has  a 
tendency  to  stop  identifying  P waves.  The  results  from  the 
computer  program  failures  show  that  more  thorough  knowledge 
of  the  relationships  between  the  P wave  and  other  wave  shapes 
on  the  ECG  need  to  be  known. 


LOCATION  OF  DISASSOCIATED 
P WAVE  IN  AN  ELECTROCARDIOGRAM 

I.  Introduction 


Background 

When  a patient  is  admitted  to  a Cardiac  Care  Unit  (CCU) , 
his  electrocardiogram  (ECG)  is  monitored  at  different  time 
intervals  depending  upon  the  severity  of  the  diagnosed  he.Vi.-t 
problem.  Until  recently,  this  monitoring  was  done  by  nurses 
in  the  CCU.  Computer  programs  which  can  do  the  monitoring 
are  now  available.  These  programs  are  used  to  inform  and 
prepare  the  nursing  staff  for  changes  toward  a worsening 
patient  condition. 

An  ECG  is  an  analog  recording  of  electrical  potentials 
from  the  heart  as  it  depolarizes  and  repolarizes.  The  im- 
pulse to  contract  begins  in  the  S-A  node  (See  Figure  1) , 
and  travels  around  the  atrium,  or  upper  chamber  of  the  heart, 
to  the  A-V  node.  The  A-V  node  is  the  only  conduction  point 
for  this  impulse  to  the  ventricles,  or  lower  chambers  of  the 
heart.  At  this  node,  the  impulse  is  delayed  and  then  con- 
ducted down  the  A-V  bundle  to  the  Purkin-je  fibers,  which 
transmit  the  impulse  to  the  remainder  of  the  heart. 

The  different  wave  shapes  on  an  ECG  correspond  to  the 
depolarization  or  repolarization  of  the  atrium  and  ventricles. 
The  QRS  complex  is  the  depolarization  of  the  ventricles  prior 
to  contraction.  The  T wave  is  the  repolarization  of  the 
ventricles.  The  wave  form  depicting  the  repolarization  of 
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the  atria  ia  hidden  by  the  QRS  complex.  In  a normal  EGG, 
as  illustrated  in  Figure  2,  the  P-8  segment  is  a relatively 
fixed  interval.  Because  the  A~V  node  is  the  only  conduction 
path  for  ho  contraction  impulse  to  the  ventricle,  the  P-R 
segment  will  lengthen  if  this  node  starts  delaying  the 
signal  more  than  normal  or  if  the  nod*  becomes  completely 
blocked.  When  this  interval,  usually  ahout  0.16  seconds, 
lengthens  to  a value  above  0.20  seconds,  the  patient  has  a 
first  degree  block  (see  Figure  3 Cal) . When  the  length  in- 
creases to  0,25  to  0.45  seconds,  the  A-V  node  is  suffici- 
ently blocked  to  completely  stop  an  occasional  impulse. 

This  will  cause  the  EGG  to  have  a P wave  with  no  following 
QRS  complex.  When  this  occurs,  the  patient  is  said  to  have 
a seccnd  degree  block  (see  Figure  3(b)).  A third  degree 
block  occurs  when  the  A-V  node  is  completely  blocked,  and  the 
atrium  and  ventricles  are  beating  at  two  different  rates. 


ECU's  showing  the  three  degrees  of  block  (a) 
first-  degree  block,  lb)  second  degree  bice’ 

(o)  third  degree  block  or  disassociated  P wave. 
(Ref  4:215) 


Fig.  3. 


This  is  termed  a disassociated  P wave  or  heart  blockage 
(see  Figure  3(c)),  depending  on  other  factors  (Ref  4:215). 

At  the  present  time  there  are  computer  programs  designed 
to  warn  of  many  cardiac  problems  (Ref  3,  6,  9) . Appendix  A 
is  a table  of  published  information  on  real-time  rhythm 
monitors.  However,  there  is  no  published  system  to  detect 
P waves  (Ref  2:293) . 

Two  previous  students  at  the  Air  Force  Institute  of 
Technology,  Flick  (Ref  1)  and  Pairett  (Ref  7) , developed 
programs  to  detect  the  P wave,  but  both  algorithms  failed 
to  find  the  disassociated  P wave.  Pairett  assumed  that  the 
P wave  would  be  directly  in  front  of  the  QRS  complex.  His 
program  achieved  96%  accuracy  in  finding  the  P wave  as  long 
as  the  P wave  was  associated.  Flick  did  not  make  the  P 
wave  association  assumption.  He  describes  two  programs  in 
his  thesis.  The  first  program  selected  a set  of  waves  for 
a feature  selection  process.  This  program  also  obtained 
frequency  domain  data  from  the  P wave,  T wave,  and  other 
sections  of  th«'  ECG.  The  second  program  located  the  P wave 
using  frequency  domain  and  rhythm  analysis.  Flick's  pro- 
gram had  a success  rate  of  86%, 

Purpose 

This  thesis  has  a two-fold  purpose.  The  first  purpose 
is  to  collect  ECG  data  from  hearts  that  lie  on  the  spectrum 
from  nearly  normal  to  cardiac  failure.  Cardiac  failure  is 
defined  in  Guyton's  Textbook  of  Medical  Physiology  (Ref  4: 
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332)  as  "failure  of  the  heart  to  pump  blood  adequately." 

The  second  purpose  of  this  thesis  is  to  develop  a computer 
program  to  find  the  disassociated  P wave.  This  program  would 
require  minimal  intervention  from  an  operator  and  would  not 
be  based  on  identification  of  the  P wave  or  its  location  with 
respect  to  other  ECG  wave  shapes. 

Assumptions 

The  ECG  for  a heart  proceeding  from  near-normal  to  death 
was  taken  from  a dog.  This  was  done  with  the  knowledge  that 
a dog's  heart  produces  an  ECG  similar  to  a human's.  It  was 
assumed  that  the  dog's  increased  heart  rate  compared  to  a 
man's,  would  not  hide  any  essential  data.  Cardiac  failure 
was  induced  by  blockage  and  an  overdose  of  drugs  on  the  as- 
sumption that  an  ECG  produced  by  these  two  techniques  would 
be  representative  of  ECG's  recorded  in  a CCU. 

The  computer  program  is  based  on  the  following  as- 
sumptions : 

1.  The  heart-beat  rate  does  not  change  faster  than 
0.01  seconds  over  a period  of  three  seconds. 

2.  This  program  is  a subroutine  for  a larger  ECG 
monitoring  program,  and  the  main  program  and  this 
program  can  be  made  compatible. 

3.  There  is  not  any  important  information  in  the  ECG 
wave  forms  above  100  Hertz. 


4 .  When  this  program  is  initiated  the  patient  is  in 
good  enough  condition  to  produce  two  consecutive 


heart  beats  without  PVC's  and  the  patient  does  not 
have  a third  degree  blockage. 


Approach 

To  gather  the  data,  two  techniques  were  used.  In  the 
first  technique  the  dog  was  anesthetized  and  a catheter  was 
inserted  in  an  artery  feeding  the  heart  muscle.  This  pro- 
cedure produced  a blockage  type  ECG.  In  the  second  technique 
the  dog  was  also  anesthetized,  but  an  overdose  of  Cedilanid-D, 
digitalis,  was  given.  This  caused  a different  type  of  ECG. 

The  ECG's  were  digitized  and  prepared  for  computer  analysis. 

The  computer  program  to  analyze  the  data  finds  the  P, 

QRS , and  T waves  in  the  first  two  heart  beats.  From  this 
identification,  the  P-P  interval,  the  distance  between  P 
waves,  is  found.  This  P-P  interval  is  then  "stepped"  on 
through  the  program.  An  interval  around  the  expected  P wave 
is  tested  for  a peak  value.  Then  another  test  is  initiated 
to  see  if  the  area  contains  a PVC  or  some  other  abnormality 
that  would  produce  an  erroneous  answer.  Finally,  a graph  is 
drawn  of  the  wave  shapes  and  the  identified  P wave. 
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XI.  Data  Collection 


General  Preparations 


There  was  a strong  need  to  have  a collection  of  electro- 
cardiograms (ECG)  proceeding  from  near  normal  heart  beat 
patterns  to  cardiac  failure.  This  collection  would  contain 
all  three  stages  of  A-V  blockage.  Premature  Ventricular  Con- 
tractions (PVC) , and  as  many  other  heart  beat  abnormalities 
as  possible.  This  data  would  serve  as  a data  library  for 
any  work  done  in  the  study  of  ECGs. 

The  ECG  library  at  Cox  Heart  Institute  had  been  searched, 
but  a suitable  collection  of  ECGs  was  not  found.  Therefore, 
the  collection  had  to  be  generated.  A dog's  heart  beat  is 
similar  to  a human's  except  it  beats  at  a faster  rate.  Hence, 
a cardiac  failure  could  be  induced  into  a dog  and  his  ECG 
recorded  on  magnetic  tape  for  future  processing,  storage, 
and  analysis.  Appendix  B shows  a block  diagram  of  the  equip- 
ment used  to  record  the  ECG. 

The  two  dogs  used  in  this  experiment  were  obtained 
through  Wright  State  Medical  School  research  channels. 

Appendix  C is  a table  showing  what  procedures  were  done  to 
each  dog  prior  to  its  arrival  at  Cox  Heart  Institute.  One 
dog  weighed  34.1  pounds  and  the  other  weighed  39.2  pounds. 

Both  were  male  German  Sheperds . 

Prior  to  the  experiment  the  dogs  were  washed  and  the 
intestinal  tract  emptied.  Each  dog  was  anesthetized  with 
Sodium  Pentothal,  and  shaved  on  the  right  hind  leg  and  left 
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side  of  the  neck,  hafore  being  brought  into  the  experiment 
room.  The  dosage  of  Sodium  Pentothal  depended  on  the  weight 
of  the  dog.  Too  much  would  kill  the  dog  and  too  little  would 
not  anesthetize  him  sufficiently.  Throughout  the  experiment, 
more  Sodium  Pentothal  was  administered  if  necessary. 

The  recording  was  done  on  an  AMPEX  portable  analog 
recorder  using  one  inch  magnetic  tape.  Only  three  channels 
were  used.  Channel  one  recorded  the  ECG  from  a chest  lead. 
Channel  two  recorded  the  electrical  signal  from  a lead  placed 
in  an  artery  beside  the  heart  muscle.  Channel  three  was  used 
as  an  editing  signal  for  data  retrieval. 

The  chest  lead  location  was  chosen  to  give  maximum  P 
wave  amplitude  and  to  conform  as  closely  as  possible  to  the 
positioning  of  ECG  leads  on  a human  in  a Cardiac  Care  Unit 
(CCU) . On  one  dog  the  maximum  P wave  amplitude  was  negative. 
This  would  serve  to  test  any  monitoring  program  l , „nst  a 
negative  P wave,  since  this  form  is  used  in  CCUs  under  cer- 
tain circumstances.  The  chest  lead  wave  shape  was  monitored 
on  a cathode  ray  tube  (CRT)  throughout  the  experiment  (see 
Figure  4) . 

The  research  team  led  by  Mr.  Karl  Kardenat  from  Wright 
State  Medical  School  at  Cox  Heart  Institute  performed  all  of 
tiie  surgery  in  this  thesis.  This  surgery  and  all  dog  prepa- 
rations used  Cox  Heart  Institute  facilities. 

The  lead  inside  the  body  was  made  by  putting  a small 
piece  of  solder  on  the  end  of  a thin  strand  of  insulated 
wire  (see  Figure  5) . This  wire  was  then  put  into  a catheter 
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Fig.  5.  The  tip  of  the  internal  heart  lead 
was  a small  piece  of  solder  on  the 
end  of  wire  in  a catheter  sleeve. 


sleeve  to  help  direct  it  down  the  artery  (see  Figure  6) . The 
a longitudinal  incision  was  made  over  the  carotid  artery  and 
the  artery  loosened  from  the  surrounding  tissue.  Next,  a 
small  cut  was  made  in  the  artery,  and  the  artery  clamped 
above  the  cut.  The  tip  of  the  lead  was  inserted  and  guided 
down  the  artery  using  a fluoroscope  x-ray  to  see  its  place- 
ment (see  Figure  7)  . The  arteries  cannot  be  seen  on  the 
fluoroscope  screen,  therefore  a technique  of  remote  directing 
was  used.  The  signal  from  this  lead  was  displayed  on  the  CRT 
monitor  (see  Figure  4).  The  lead's  placement  was  adjusted  to 


E 


Fig.  6.  The  complete  internal  lead  and  catheter  system  is  shown. 


Fig.  7,  The  placement  of  the  internal  lead  was  monitored  using 

a fluoroscope  x-ray.  The  faint  thin  line  is  the  internal 
lead.  The  round  spot  with  a dark  line  coming  from  it 
is  the  external  chest  lead. 
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give  a maximum  P wave  spike.  Once  this  lead  was  in  place,  the 
artery  was  clamped  around  the  lead  wire  to  minimize  blood 
loss.  The  wave  form  from  this  lead  gave  an  amplified  P wave 
compared  to  that  seen  on  an  ECG.  In  one  case  the  wave  form 
corresponding  to  the  P wave  was  higher  than  the  remaining 
wave  shapes  on  the  ECG.  This  lead  would  help  point  out  the 
P waves  on  the  ECG. 

The  third  channel  was  an  editing  channel  where  a one  volt 
signal  on  the  tape  would  indicate  a wave  form  to  be  digitized 
in  the  other  data  channels.  This  voltage  came  from  a direct 
current  power  supply  through  a double-pole  single  throw 
switch.  The  input  was  monitored  by  a digital  voltmeter  at 
the  recorder  input. 

Also  displayed  on  the  CRT  was  the  systemic  blood 
pressure  (SEP)  (see  Figure  4} . This  wave  form  was  collected 
by  a transducer  actuated  by  fluid  pressure  from  the  blood  in 
an  artery  in  the  right  hind  leg.  This  wave  shape  was  very 
important  in  monitoring  the  condition  of  the  dog  during  the 
experiment. 

With  all  the  preparation  done,  the  dog  was  ready  for  the 
experiment.  Two  means  of  inducing  cardiac . failure  were  used: 
one  on  each  of  the  two  dogs.  The  first  means  was  through 
blockage  of  the  artery  feeding  the  A-V  node  section  of  the 
heart  muscle.  The  second  method  was  through  an  overdose  of 
a toxic  drug. 
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Elockage  Technique 

The  blockage  technique  used  a catheter  with  a coil  of 
copper  wire  placed  at  the  end  of  the  inner  catheter.  The 
entire  catheter  system  is  made  up  of  an  outer  radiopaque 
(Formacath)  polyethylene  catheter  with  an  outer  diameter  of 
approximately  Q.11Q  inches  and  a length  of  about  11.8  inches. 
The  vinyl  inner  catheter  had  an  outer  diameter  of  approxi- 
mately 0.Q65  inches  and  was  19.7  inches  long.  The  back  end 
of  the  inner  catheter  had  a three-way  stopcock  attached. 

The  outer  catheter  had  an  adapter  connected  to  its  back  end 
that  would  allow  the  inner  catheter  to  move  with  a slight 
amount  of  friction.  The  front  end  of  the  catheter  was  bent 
at  a right  angle  to  the  coronary  ostium  plane,  after  the 
catheter  was  put  at  the  root  of  the  aorta.  The  inner  catheter 
was  then  placed  through  the  lumen  of  the  outer  catheter  (see 
Figure  8) . 

The  blockage  device  was  a solid  copper  wire  0.02  inches 

in  diameter  that  was  bent  to  form  a helical  shape.  The 

length  of  the  helical  wire  was  approximately  Q.Q4  inches. 

The  length  of  the  coil  helped  determine  the  clotting  time 

for  a formation  around  the  wire  and,  thus,  the  time  to  cardiac 

\ 

failure.  From  experiments  performed  by  R.  Karl  Kordenat,  this 
time  could  vary  from  an  hour  to  several  days  (Ref  5:363)  . 

The  coil  is  placed  on  the  front  end  of  the  Inner  catheter. 

The  coil  should  have  some  resistance  from  slipping  off  but 
still  be  loose  enough  for  the  outer  catheter  to  push  it  off. 


Once  the  catheter  system  is  prepared,  it  is  placed  in 
the  same  opening  in  the  carotid  artery  as  the  internal  ECG 
lead  wire.  The  elasticity  of  the  artery  wall  is  sufficient 
and  the  ECG  wire  small  enough  to  make  this  possible.  Once 
again,  under  fluoroscopic  control,  the  catheter  is  snaked 
through  the  artery  system  to  the  tip  of  the  coronary  ostium. 
The  catheter  is  then  directed  into  the  anterior  descending 
branch  of  the  coronary  ostium.  A contrast  media  (Hypoque  75) 
is  periodically  injected  through  the  inner  catheter  to  help 
show  the  artery  locations. 

Once  the  tip  of  the  inner  catheter  is  in  the  desired 
position,  the  outer  catheter  is  pushed  forward,  forcing  the 
coil  off  the  inner  catheter.  The  coil  is  now  free  in  the 
artery  and  can  he  "directed"  to  the  desired  position.  Like 
the  size  of  the  wire  coil,  the  position  of  the  coil  in  the 
artery  with  respect  to  the  branches  off  the  artery  is  impor- 
tant for  determining  the  time  until  cardiac  failure.  The 
position  is  also  important  in  determining  what  wave  forms 
the  heart  will  produce  prior  to  failure  (Ref  5:363) . For 
example,  if  placed  too  high  in  the  anterior  descending  branch 
P wave  disassociation  has  a reduced  chance  of  happening. 

Once  the  coil  was  in  place,  the  heart  gave  good  data 
for  55  minutes  before  sinus  rhythm  was  lost  and  could  not  be 
regained.  PVC's  were  observed  about  3Q  minutes  after  the  coil 
was  released.  Xylocaine  and  Pronestyl  were  administered  to 
help  the  heart  recover.  A few  minutes  later  Sodium  Bicar- 
bonate was  also  given  to  help.  These  drugs  reduced  the 


occurrence  of  tiie  PVC's  and  the  heart  began  beating  more  in 
sinus  rhythm.  About  13  minutes  later  the  PVC’s  returned  and 
drugs  had  no  effect  on  them. 

Toxicity  Technique 

For  this  procedure  the  dog  was  prepared  as  described 
above.  However,  instead  of  inserting  a helical  coil,  Cedi- 
lanid-D,  digitalis,  was  injected  into  the  bloodstream.  The 
weight  of  the  dog  determined  the  dosage.  An  injection  of 
Q.8  mg  was  made  every  ten  minutes  until  a depressed  S-T 
segment  on  the  ECG  was  observed.  This  was  observed  15  minutes 
after  the  first  injection.  Following  the  initial  depressed 
S-T  segment,  three  more  injections  were  made  at  approximately 
one  hour  after  initial  injection.  No  drugs  were  given  to 
prolong  the  heart  action  because  of  their  reaction  with  the 
digitalis. 

Analog  to  Digital  Conversion 

The  data  on  the  analog  tape  had  to  be  put  into  a format 
which  the  digital  computer  would  accept.  To  accomplish  this, 
the  analog  tape  was  played  back  on  a Bell  and  Howell  VR370QB 
tape-deck  set  on  the  intermediate  band  width  at  a speed  of 
7-1/2  inches  per  second.  This  output  was  sent  to  a Comcore 
5Q0Q  analog  computer.  This  computer  multiplied  the  incoming 
signal  amplitude  by  50  and  then  passed  the  signal  through  a 
capacitor  to  remove  the  base  line  drift.  The  signal  was  then 
passed  through  a 100  Hertz  (Hz)  low  pass  filter.  This  filter 
reduce  1 signals  above  1Q0  Hz  at  a rate  of  3 dB  per  octave. 
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This  step  waa  designed  to  remove  moat  of  the  noise  in  the 
signal.  At  this  point  the  signal  path  was  split. 

One  signal  went  to  a strip  chart  that  gave  a continuous 
readout  of  the  signal.  The  second  signal  went  to  the  Xerox 
Sigma  7 computer.  This  computer  digitized  the  signal.  The 
computer  sampled  the  analog  data  every  two  milliseconds  using 
16  bits  per  sample.  The  bit  precision  exceeded  the  recom- 
mended digitization  rate  of  500  samples  per  second  using  9 
bits  per  sample  set  by  the  Committee  of  Electrocardiography 
of  the  American  Heart  Association  (Ref  8) . The  number  of 
data  points  per  block  of  data  (512)  was  chosen  to  fill  the 
Control  Data  Corporation,  CDC  6600  computer's  buffer. 

The  digitized  signal  was  then  recorded  on  an  on-line 
digital  tape  system.  The  digitized  signal  was  also  routed 
back  to  the  strip  chart  through  a digital  to 'analog  converter. 
This  produced  a read-out  on  the  strip  chart  of  both  the  con- 
tinuous data  and  the  sampled  data. 

An  existing  program  was  used  to  convert  the  data  from 
Sigma  7 'a  32  bit  word  to  CDC  6600's  16  bit  word.  The  output 
of  this  conversion  was  written  onto  a magnetic  tape.  A pro- 
gram was  then  written  to  change  the  data  on  this  tape  to  an 
eight  digit  field  with  an  accuracy  of  two  decimal  places. 

This  program  is  in  Appendix  C. 
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III.  Computer  Program 


Background 


The  previous  approaches  used  to  solve  the  P wave  iden- 


tification problem  used  a function  that  transformed  the  wave 


shape  into  a different  domain  (Ref  1,  2).  After  comparing 


the  wave  shape  in  Figure  9(a)  to  that  in  Figure  9(b),  there 


seems  to  be  no  resemblance  between  them.  Therefore,  com- 


paring the  different  domain  characteristics  would  lead  to 


false  identification. 


r r 
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Fig.  9.  This  ECG  is  from  a heart  progressing 
from  near  normal  (a)  to  disassociated 
P wave  (b) . 


One  important  characteristic  of  the  P wave  is  its 


periodicity.  The  distance  between  P waves,  P-P  interval,  at 


P wave  di.sassoeiation  had  changed  0.07  second  from  the  P-P 


interval  in  a normal  heart.  The  cardiologist  typically  uses 


this  characteristic  to  identify  the  P waves.  He  finds  the 


R wave  and  then  looks  for  a peak  above  the  noise  about  0.16 


seconds  before  the  R wave  and  calls  this  peak  a P wave. 

This  procedure  is  repeated  for  the  next  heart  beat.  One 
point  of  a pair  of  dividers  is  placed  on  the  peak  of  each 
P wave,  giving  the  P-P  interval. 

With  this  interval  known,  the  physician  moves  the 
dividers  down  the  ECG  looking  for  wave  form  peaks  around 
each  point.  If  the  peaks  become  consistently  different  from 
where  the  divider  point  is,  the  heart  rate  has  become  faster 
(smaller  P-P  interval) ; or  become  slower  (longer  P-P  inter- 
val) . Then  the  interval  represented  by  the  dividers  must  be 
adjusted. 

When  the  ECG  contains  a disassociated  P wave,  the  inter- 
val on  the  dividers  is  virtually  the  only  way  to  identify  the 
P wave.  The  difficulty  is  that  during  disassociation,  there 
are  usually  many  other  wave  forms  present;  such  as,  Premature 
Ventricular  Contractions  CPVC) , which  cover  or,  at  best,  make 
the  P wave  less  apparent  (see  Figure  9 (b) ) . Frequently  the 
P wave  will  not  be  apparent  for  several  beats  because  it  is 
"covered  up"  by  wave  shapes  having  higher  amplitudes.  In 
this  case  the  dividers  are  just  moved  down  the  ECG  as  normal 
until  a P-P  interval  can  be  identified.  Then  the  distance  on 
the  dividers  is  adjusted,  and  this  distance  is  moved  back  up 
the  ECG  looking  for  P waves. 

This  technique  depends  heavily  on  having  at  least  two 
associated  P waves  before  disassociation  occurs,  or  the  task 
becomes  significantly  harder.  Even  when  having  these  two  P 
waves  present,  a long  period  of  time  can  be  required  to 


find  these  P waves.  This  time  requirement  complicates  the 


real-time  analysis  of  ECG's. 


Initializing  the  Program 


This  program  is  patterned  after  the  approach  used  above. 


The  program's  list  of  variables,  a copy  of  the  program,  and 


explanation  of  the  program  by  sections  is  in  Appendix  E.  The 


main  points  and  unique  elements  of  the  program  will  now  be 


covered . 


The  program  is  designed  to  require  a minimum  amount  of 


interaction  with  the  operator.  However,  the  operator  does 


have  to  enter  the  number  of  peaks  in  the  ECG.  This  will  come 


from  looking  at  the  patient's  ECG  strip.  If  the  patient's 


heart  beat  has  joined  the  P wave  of  the  next  beat,  there  will 


be  three  peaks.  If  the  T and  P have  joined,  there  will  be 


two  peaks.  The  operator  will  also  need  to  enter  whether  the 


P wave  is  positive  or  negative.  This  information  comes  from 


the  patient's  ECG.  Whether  the  P wave  is  positive  or  nega- 


tive, depends  on  which  ECG  lead  placement  scheme  is  used 


(Ref  4:195).  If  the  P wave  is  positive,  a zero  (0)  is 


entered;  if  the  P wave  is  negative,  a one  (1)  is  entered. 


The  operator  will  also  enter  which  block  of  data  he  first  wants 


to  see,  and  how  often  he  wants  to  see  the  output. 


Reading  the  Data 


Next  the  program  calls  the  subroutine  READ  to  read  the 


digitized  wave  shape  values  in  blocks  of  512  data  points. 


After  the  data  is  read  into  memory,  a matrix  of  integer 
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data  values  is  formed.  This  matrix  is  formed  by  ignoring  all 
the  values  below  a whole  number.  A rounding  off  process  is 
not  used  because  the  exact  maximum  amplitude  value  is  not 
important.  If  the  P waves  are  positive,  all  negative  data 
points  are  s«'t  equal  to  zero..  If  the  P waves  are  negative, 
the  absolute  value  of  each  data  point  is  found.  Then  a 
threshold  value  below  the  peak  value  of  the  P wave  is  found, 
and  all  data  below  this  threshold  is  set  equal  to  zero.  These 
two  processes  are  designed  to  eliminate  any  noise  spikes 
which  would  be  identified  as  maximum  values  in  the  first 
program  task. 

Initial  P-P  Interval 

After  the  data  has  been  initially  processed,  a decision 
is  made  whether  this  is  the  first  data  block  or  not.  If  this 
is  not  the  first  data  block,  a P-P  interval  has  been  found 
and  control  is  transferred  further  down  the  program  to  pro- 
cess the  data.  If  this  is  the  first  data  block,  the  first 
two  P waves  must  be  found  and  then  the  P-P  interval  deter- 
mined. 

To  find  the  first  two  P waves,  a subroutine  MAX  is 
called  that  identifies  the  maximum  values  in  the  first  two  heart 
beats.  The  basic  idea  is  that  the  slope  of  the  line,  amplitude 
difference  between  two  adjacent  points,  is  found  and  compared 
to  the  last  slope.  The  three  conditions  are  possible:  slopes 
are  positive  (upward);  negative  (downward);  and  zero  (flat). 

All  changes  in  slope  are  labeled  as  follows;  a change  to  a 
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positive  slope  is  labeled  one  Cl 1 ; a change  to  a negative 
slope  is  labeled  zero  (Q) ; and  a zero  slope  does  not  change 
the  slope  label.  The  maximum  point  is  the  point  before  the 
slope  label  changes  from  one  to  zero. 

Figure  IQ  is  an  example  of  how  this  subroutine  works. 

A value  of  one  Cl)  is  initially  assigned  to  the  last  slope 
variable,  LSIGN . The  amplitude  difference  between  the  last 
data  point  and  the  present  data  point  is  then  found.  If 
this  difference  is  positive,  the  present  slope  variable, 
JSIGN,  is  assigned  a one.  The  process  continues  with  the 
remaining  data  points  as  long  as  the  slope  is  positive.  If 
the  slope  changes  to  zero,  as  in  Figure  10  Cb)  and  (c) , the 
values  of  LSIGN  and  JSIGN  remain  the  same.  When  the  slope 
changes  to  negative  (see  Figure  10(d),  JSIGN  changes  to  zero. 
Since  this  is  different  from  LSIGN  (which  is  1) , a maximum 
point  is  defined  as  the  last  point  where  JSIGN  and  LSIGN 
equalled  one  (1)  , (MAX  1 in  Figure  10).  The  subroutine  con- 
tinues until  all  the  data  in  the  data  search  window,  as  de- 
fined in  other  parts  of  program,  has  been  searched  for  maxi- 
mum points.  If  more  than  one  maximum  point  is  found  (see 
Figure  10) , a check  is  initiated  to  see  how  close  the  maxi- 
mum points  are.  If  they  are  within  a set  interval  (3  data 
points  in  this  program) , the  maximum  value  of  these  points 
is  found  and  called  the  maximum  point. 

After  the  maximum  points  are  found,  the  main  program 
identifies  the  R wave  in  the  first  heart  beat  by  finding  the 
maximum  point  of  all  the  maximum  values  returned. 
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The  maximum  point  before  the  R wave  peak  is  identified  as  the 
P wave.  If  the  R wave  is  the  first  wave  form,  the  last  wave 
form  in  the  heart  beat  is  called  the  P wave.  Once  the  wave 
forms  in  the  first  heart  beat  are  found,  the  process  is  re- 
peated for  the  second  heart  beat. 

Once  both  P waves  are  identified,  the  distance  between 
them  is  the  number  of  data  points  found  and  is  called  the 
P-P  interval.  Since  the  sample  rate  is  500  samples  per 
second,  this  P-P  interval  could  be  multiplied  by  .002  to 
find  the  interval  in  time. 

Remaining  P-P  Intervals 

With  the  initial  P-P  interval  found,  the  subroutine 
GREATER  is  used  to  find  the  remaining  P waves.  Subroutine 
GREATER  defines  certain  parameters  for  subroutine  MAX  to  use 
to  find  the  next  P wave.  These  parameters  define  where  to 
start  looking  for  the  next  P wave  and  the  data  search  window. 
Initially,  this  window  is  set  to  ±2  data  points  from  the 
expected  P wave  location.  If  a peak  is  not  found,  the  data 
window  is  expanded  by  two  (2)  data  points  in  both  directions, 
until  the  window  reaches  a maximum  value  of  30  data  points. 
The  search  for  peak  values  is  accomplished  by  subroutine 
MAX.  If  more  than  one  maximum  value  is  returned,  the 
largest  value  is  used.  If  no  peak  is  found,  the  P wave  is 
assumed  to  be  where  it  was  expected.  This  assumes  the  heart 
rate  will  not  change  significantly  until  a P wave  is  found. 
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The  location  of  the  P wave,  or  where  it  should  have  been  (if 
a peak  was  not  found) , is  transferred  to  the  main  program. 

The  value  transferred  from  subroutine  GREATER  is  used 
to  find  the  new  P-P  interval.  If  the  P wave  was  not  at  the 
expected  location,  a different  P-P  interval  is  found  and 
placed  in  the  P-P  interval  matrix.  If  the  P wave  was  at  the 
expected  location  or  peak  was  not  found,  the  P-P  interval 
does  not  change  from  the  previous  one;  and  that  value  is 
placed  into  the  P-P  interval  matrix.  Figure  11  is  an 
example  of  the  interval  remaining  the  same  and  then  changing. 

Output 

The  output  of  the  program  is  accomplished  by  a call  to 
subroutine  MARK,  which  in  turns  calls  subroutine  AXIS  and 
PMARK.  Subroutines  were  used  to  facilitate  changing  methods 
of  output.  This  particular  program  uses  a Textronix  4014 
terminal  for  CRT  graphing  with  a Textronix  4631  copier  at- 
tached to  the  terminal  for  paper  copies  of  the  CRT  display. 

Subroutine  MARK  initializes  or  erases  the  screen;  calls 
AXIS;  plots  a data  block;  sets  parameters  for  PMARK;  and 
prints  the  value  of  the  P-P  intervals  in  that  data  block. 

AXIS  is  designed  to  draw  a frame  around  the  output  presented. 
Then  the  X and  Y axes  are  drawn  and  divided  by  tic  marks . 

This  subroutine  also  labels  the  major  tic  marks  and  both  axes. 
After  all  of  this  is  drawn,  control  is  returned  to  sub- 
routine MARK. 
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hanged  interval . 


MARK  now  plots  the  data  and  checks  to  see  if  a plotted 
data  point  corresponds  to  a peak  of  an  identified  P wave. 

If  the  data  point  does  correspond,  parameters  for  PMARK  are 
set,  and  PMARK  is  called. 

PMARK  draws  a straight  vertical  line  from  the  identified 
P wave  peak.  This  line  is  used  to  easily  identify  the  wave 
forms  selected  by  the  program.  Once  the  P waves  have  been 
marked,  control  is  again  returned  to  subroutine  MARK,  which 
now  prints  out  the  P-P  intervals  in  the  block  of  data  being 
shown  on  the  screen.  After  this,  the  program  is  halted  so 
that  the  graph  can  be  studied  or  copied.  When  the  program 
is  started  again,  the  screen  is  prepared  for  the  next  plot. 

Continuing  the  Main  Program 

After  the  output  has  been  generated,  control  is  re- 
turned to  the  main  program.  Now  the  operator  is  given  an 
opportunity  to  exit  the  program,  or  have  another  block  of 
data  read  into  memory.  If  the  operator  elects  to  continue 
the  program,  another  block  of  data  is  read.  A system  to 
transfer  the  area  of  the  next  expected  P wave  is  needed  to 
be  devised.  Since  there  are  almost  1.5  million  data  points 
to  be  processed,  the  data  has  to  be  analyzed  in  block  seg- 
ments. This  requires  a system  to  carry  forward  into  the 
next  data  block,  a part  of  the  P-P  interval  not  used  in  the 
previous  segment.  For  example  (refer  to  Figure  11  for 
assistance),  if  the  P*-P  interval  is  291  data  points,  but 
only  10Q  data  points  are  left  in  this  data  block,  then  the 
first  191  data  points  must  be  ignored  in  the  next  data  block 
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to  get  the  data  window  around  the  calculated  P wave.  This 
skipping  routine  consists  of  finding  the  number  of  data  points 
left  in  the  first  data  block  (LB  = 512  - maximum  point) , and 
then  subtracting  this  from  the  known  P-P  interval.  Once  the 
skipping  routine  is  finished,  the  next  problem  is  how  to 
calculate  the  next  P-P  interval.  This  interval  is  calculated 
by  adding  LB  to  the  number  of  data  points  corresponding  to 
the  identified  P wave.  Since  this  last  addition  routine  can 
only  be  performed  when  going  from  one  data  block  to  another, 
a decision  step  must  be  included  as  to  whether  to  use  this 
step  or  not.  The  last  addition  routine  is  used  only  when 
the  difference  between  the  P wave  data  point  is  less  than  the 
last  calculated  P-P  interval.  Therefore,  the  decision  step 
is  based  on  whether  the  location  value  is  less  than  the  ex- 
pected P-P  interval . 
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IV.  Results 


Data  Col loot  ion 

Data  wore  collected  from  three  dogs.  The  blockage  tech- 
nique was  used  on  two  dogs:  Dog  ] and  Dog  2.  The  toxicity 
technique  was  used  on  Dog  3.  The  electrocardiograms  (ECC) 
from  these  dogs  wore  recorded  on  separate  reels  of  analog 
tape  for  future  use. 

The  ECG's  from  Dog  2 and  Dog  3 were  analyzed  by  a 
cardiologist  and  disassociated  P waves  identified.  The  FCG 
from  Dog  2 had  a negative  P wave  because  the  AVL  load  was 
used  as  the  exterior  lead.  This  provided  a data  base  to 
teat  the  computer  program  against  negative  P waves.  Figure 
13(a)  is  a sample  of  this  ECC. 

The  ECC  from  Dog  3 had  a positive  P wave  (see  Figure 
13(b)),  but  there  were  abnormal  wave  shapes  at  1.6  Hz 
throughout  the  ECC.  These  wave  shapes  were  apparently  caused 
by  the  Frequency  Modulated  (FM)  analog  recorder  being  over 
modulated  by  a strong  signal  going  into  it.  This  accidental 
loss  of  signal  proved  to  be  valuable  in  order  to  test:  the 
computer  program  against  abnormal  wave  forms  in  an  KCC . A 
patient  in  a Cardiac  Care  unit  frequently  produces  artificial 
wave  forma  different  from  a normal  heart.  A computer  program 
to  find  a p wave  must  take  this  into  account.  Figure  13(c) 
is  an  example  of  Dog  3's  ECC  with  these  wave  forms. 

The  wave  forms  on  an  ECC  of  a heart  dying  trom  blockage 
is  significantly  different  from  that  of  a heart  dying  from 
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Fig.  13.  (a)  Dog  2’s  F.CG  has  a negative  P wave. 

(b)  Dog  3's  ECG  has  a positive  P wave, 

(cj  Dog  3*s  ECG  has  a data  drop  out 

occurring  at  a 1.6  Hz  rate. 


toxicity.  The  ECG  from  a blocked  heart  has  wave  forms  with 
a higher  amplitude  peaks  and  longer  duration.  The  ECG  from 
the  toxic  heart  has  lower  peaks  and  is  more  rhythmic.  In 
general,  a heart  dying  from  a blockage  produces  a more 
"violent"  ECG.  This  ECG  was  harder  for  the  cardiologist  to 
read . 

To  read  this  ECG  the  cardiologist  used  the  internal  lead 
to  help  find  the  P wave.  Generally,  this  internal  lead 
provided  a good  wave  shape  to  help  find  the  P wave.  On 
occasion,  near  the  end  of  the  ECG  and  the  death  of  the  dog, 
the  P wave,  evt  n in  this  lead,  could  not  be  found. 

The  violence  around  the  P wave  disasso.ciation  area  on 
an  ECG  is  shown  in  Figure  14 . 
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Fig.  14.  The  violence  around  the  P wave 
disassociation  experienced  on 
Dog  3 * s ECG . 


Computer  Program 

The  computer  program  uses  the  last  P-P  interval  to  find 
the  area  of  the  next  P wave.  This  demands  that  the  data  be 
handled  sequentially  without  any  breaks,  once  the  P-P  inter- 
val is  initially  calculated.  The  computer  program  processed 
the  data  with  an  output  every  50th  data  block  after  the  first 
block  of  data.  The  program  showed  a promising  ability  to 
pick  out  P waves  in  many  instances  (see  Figures  15  and  16) 
and  correct  itself  at  other  times  (see  Figure  17) . 

The  computer  program  was  stopped  when  the  first  wrong 
wave  form  was  detected.  To  find  why  the  program  failed,  the 
program  was  restarted  with  an  output  of  every  10th  data 
block  and  finally  every  data  block  until  the  cause  was  iso- 
lated. This  procedure  was  accomplished  twice  during  the 
first  50  data  blocks,  and  again  during  the  second  50  data 
blocks. 

The  first  problem  was  in  data  block  43  (see  Figure  18) . 
The  program  picked  out  the  first  P wave  in  the  block,  and 


' ft*  •>!>. • lA V"  t i 


jar.; . 

I?. , , 

fo'  ...  'Vi 


m.  • . * 


gS*.  , . 

S'. 

pVUI 


r.i  ■ 


,'li 

8?'  *.  .)l 

,hl 
* Sj 

. ] 


pi:,  l-  -ft- 


,, ' 


. . */.M® 


fv- 


0! 

• 

> -P 

ITS 

a 

S 

0 

a 

0 

as 

p 

x:  t3 

s', 


’O 

-U 

4J  <D 

•to 

3 -P 

■M 

O ns 

I 

■H 

,WY 

T3  U 

. & 

<d  o 

X w 

W 

0 M 

fjj 

•H  (fl 

a 

m 

a) 

w 

e 

m -p 

ffl 

p a 

fj| 

Di  CO 

s 

0 

p as 

1 

ax 

■p 

;|j| 

p 

sj 

as  as 

jB 

4J  P 

'« 

3 0 

u5| 

am 

l>s|| 

e a) 

tgH 

0 Si 

0 

F&9 

*•  • <•  _s, 


it  “i'W  ‘-rVT!T*V.w>. 


then  the  data  produced  wave  shapes  that  the  program  could  not 
handle.  After  that  data,  the  program  became  an  R wave  detec- 
tor. The  reason  for  the  strange  data  wave  shapes  is  unknown 
but  may  reflect  a singular  error  in  the  data  recording 
scheme . 

The  second  problem  was  in  data  blocks  74  and  75. 

Through  data  block  74  (the  graph  up  to  the  broken  line  in 
Figure  19),  the  program  picked  out  the  correct  P waves.  In 
data  block  75  the  first  P wave  was  identified  but  only  part 
of  the  second  P wave  was  present  because  of  the  1.6  Hz  data 
loss  problem  described  previously.  The  program  identified 
this  small  amount  of  data  as  the  P wave.  When  the  P-P  in- 
terval was  transferred  to  the  second  wave  shape,  the  peak 
of  the  P wave  was  outside  of  the  expanded  data  search  window. 
Therefore,  the  program  assumed  the  P wave  was  not  there  and 
the  last  P-P  interval  was  retained.  With  the  next  data 
block,  the  program  became  a T wave  detector. 


is  the  combination  of  data  blocks  74  and  75 
block  74  is  presented  up  to  the  broken  line 


V.  Conclusion  and  Recommendation 

Conclusion 

Two  main  conclusions  concern  the  precision  of  the  data 
and  the  success  of  the  computer  program.  First,  the  sample 
rate,  10QQ  samples  per  second,  was  too  high.  This  produced 
redundant  data.  The  data  rate  of  500  samples  per  second  was 
more  manageable  and  met  the  needs  of  the  program.  The 
analyzed  data  also  did  not  need  to  be  sampled  to  a 0.Q1 
millivolt  accuracy.  The  program  is  designed  to  find  the 
P wave,  not  its  absolute  peak  value.  A relative  relationship 
between  wave  form  values  was  sufficient.  The  0.01  millivolt 
accuracy  caused  problems  by  inducing  unwanted  noise  into  the 
data  which  had  to  be  taken  out.  A flat  peak  was  caused  by 
the  ignoring  of  fractional  millivolt  variations  and  the 
sampling  routine.  However,  since  the  graphing  pen,  or  CRT 
beam,  has  a finite  width,  an  absolute  peak  is  not  desired, 
this  flat  peak  is  acceptable. 

The  computer  program  demonstrated  that  an  adjustable 
data  search  window  is  needed.  When  the  search  window  was 
set  toe  large  (3Q  data  points),  whole  heart  beats  were 
skipped,  particularly  in  the  data  close  to  disassociation. 
This  was  determined  by  human  analysis  of  the  data  using  the 
computer  program  logic,  given  that  the  first  two  P waves  were 
identified  correctly.  The  problem  in  this  area  of  the  ECG 
seems  to  be  the  increased  heart  rate.  This  compresses  the 
wave  shapes,  putting  the  T and  R waves  within  the  large 
search  window.  With  these  wave  shapes  inside  the  window, 
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the  program  picked  the  highest  peak  (the  R wave)  as  the  P 
wave.  With,  the  data  window  small  (5  data  points),  the  P wave 
was  easily  missed  during  the  early  part  of  the  data.  This 
was  because  the  wave  forms  were  not  compressed  and  large. 

The  0.07  second  change  in  the  P-P  interval  cited  previously 
represents  35  data  points,  which  is  a significant  change  when 
wave  forms  become  close  to  each  other.  The  solution  was  an 
adjustable  data  search  window  that  started  at  5 data  points 
and  expanded  to  a maximum  of  30  data  points  if  a peak  was  not 
previously  identified. 

The  computer  program  failures  could  not  be  corrected 
because  only  one  piece  of  information,  periodicity  of  the  P 
wave,  was  used  to  identify  the  P wave.  Many  other  factors 
such  as  present  wave  shape,  relationship  to  the  R wave,  and 
historical  shape  of  the  P wave  are  used  by  a cardiologist  to 
identify  the  P wave.  These  factors  will  be  especially  im- 
portant in  identifying  disassociated  P waves . While  this 
program  failed  to  identify  P waves  at  human  accuracy  rates, 
it  developed  another  useful  approach  that  can  be  added  to 
existing  techniques  to  find  the  P wave. 


Recommendation 

The  first  thing  to  be  studied  is  the  outline  of  a thorough 
set  of  relationships  between  the  P wave  and  the  other  wave 
shapes  on  an  ECG.  At  the  present  these  relationships  are  not 
well  defined.  The  data  produced  by  this  thesis  can  be  useful 
in  defining  these  relationships.  However,  care  must  be 
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exercised  in  applying  findings  based  on  this  data  because 
there  are  small,  but  important-,  differences  between  the  dog 
and  the  human  ECG. 

After  these  relationships  are  defined,  they  should  be 
thoroughly  tested  against  human  data.  This  will  probably  re- 
quire certain  constants  and  some  rules  to  be  adjusted.  Once 
these  rules  are  thoroughly  tested,  a computer  program  can  be 
written  to  implement  the  rules. 

If  this  program  produces  a sufficiently  high  probability 
of  being  correct,  a digital  hardware  package  could  be  de- 
signed to  implement  the  program.  This  hardware  should  be  in 
a small  package  using  the  latest  in  solid  state  technology 
to  analyze  ECG’s  on  a real  time  basis.  The  small  size  would 
facilitate  each  patient  in  a CCU  having  his  own  monitor, 
without  adding  appreciably  to  the  equipment  around  his  bed. 
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Fig.  20.  Block  diagram  of  the  equipment 
used  to  record  ECGs  . 
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Table  2.  Data  on  Dogs 


Dog  1 Dog  2 


Weight 

Breed 

Animal  Procedures: 

Dipped  for  Ectoparisites 

DHL.  Vaccine 

Rabies  Vaccine 

Voice  Biopsy 

Treated  for  Ear  Mites 

Microfilaria  Exam 


34.1  lbs. 

39.2  lbs 

German 

Shepard 

German 

Shepard 

Yes 

Yes 

Yes 

Yes 

Yes 

Yes 

Yes 

Yes 

Yes 

Yes 

Negative 

Negative 
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PROGRAM  COMP  (INPUT, OUTPUT, TAPE?, TAPE21 
DIMENSION  DATA  (5121 
DO  3 1=1,2786 
READ (7)  (DATA (J) ,J=1, 5121 
IF(EOF(7)  .NE.Ol  GO  TO  900 
PRINT (2,131  (DATA(K) ,K=i,512) 
13F0RMATC10F8.21 
3C0NTINUE 
90QSTOP 
END 


Fig.  21.  Program  used  to  change  data  format. 


Figure  21  shows  the  computer  program  used  to  change  the 
data  into  a format  usable  in  the  P wave  identification  pro- 
gram after  the  original  data  was  put  into  CDC  6600  word 
length.  The  program  reads  the  data  in  an  unformated  form, 
changes  the  format  to  floating  point  with  a field  length  of 
eight  characters  including  two  decimal  places.  This  new 
format  is  then  written  onto  a drum  storage  unit  for  storage. 
This  program  can  be  expanded  to  punch  the  data  onto  cards, 
print  the  data  onto  paper,  or  use  magnetic  ‘ape  for  storage. 


Job  card. 

FTN . 

REQUEST,  TAPE2,f*PF. 

LABEL, TAPE7 , R, L=lfnl , VSN=vsn , NORING . 

LGO. 

CATALOG,  TAPE 2 , pfnl,RP=999. 

EXIT,S , 

REQUEST,  TAPE2 , *Q . 

ROUTE , TAPE2 , DC=PR, TID=id, FID=1 fn2 , ST=sv s . 


Fig.  22.  List  of  control  cards  used  to  execute 
the  program  listed  in  Figure  21. 
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Figure  22  shows  the  control  cards  necessary  to  compile 
the  program,  mount  the  magnetic  tape  containing  the  data 
converted  to  CDC  660Q  format,  request  permanent  file  space 
for  the  changed  data,  execute  the  program,  and  catalog  the 
changed  data.  A series  of  control  cards  (starting  with 
EXIT,S.)  is  also  listed  that  routes  the  output  to  the  user's 
intercom  terminal  if  the  program  should  be  terminated  before 
it  is  finished.  The  variables  used  in  these  control  cards  are 
as  follow: 

lfnl  - The  file  name  written  of  the  magnetic  tape. 

vsn  - The  number  assigned  to  the  tape  to  be  read. 

pfnl  - The  permanent  file  name  assigned  to  the  new 
data. 

id  - The  intercom  identification  letters  assigned 
to  the  operator. 

Ifn2  - The  local  file  name  assigned  to  the  routed 
file. 

sys  - The  system  designator  receiving  the  new  file. 
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In  this  appendix  the  program  to  find  P waves  will  be 
presented  and  discussed  in  detail.  The  first  part  will  be 
a detailed  discussion  of  the  program  in  short  sections. 

The  final  part  will  list  the  entire  program. 

Program  Discussion 

The  computer  program  will  be  presented  in  sections  as 
listed  in  the  controlling  program.  Since  most  calculations 
are  accomplished  in  the  subroutines,  these  subroutines  will 
be  discussed  when  they  are  first  called.  Subsequent  calls 
to  these  subroutines  will  reference  the  figure  number. 

Setting  Parameters . — The  first  few  lines  of  the  program 
(see  Figure  22)  identify  the  program,  establish  storage 


PROGRAM  FIND(TAPE2, INPUT, OUTPUT, PLOT) 

C THIS  IS  MAIN  PROGRAM  TO  FIND  3 WAVES  IN  AN  EG 5 

DATA (512) »Ni»LK,IP( 3000) ,N/ GREAT /TIP,N0/MAR</I3AUD, J 
) 

MAX/ISTEP,IMAX(f  ) ,IPP(6)  ,KK,NB,NC,M 
PP 


COMMON  I 
1, PP (273 i 
COMMON  / 
INTEGER 
REWIND  2 


c 

WRITE* » 3 
READ*,  I 

WRITE*,  2 
READ* , N 
WRITE*, 3 
IIS  NEGAT 
READ*,  N1 
WRITE* , 3 
READ* » M A 
WRITE* , 4 
READ*, M3 


READ  IN  THE  IPAUO  VALUE  FDR  I NITT 
TAKE  THE  LISTEO  BAUD  RATE  AND  DIVIDE  8Y  ID 
2HTHE  BAUD  RATE  DIVIDED  BY  13  IS 
AUO 

SET  UP  READ  AND  ANALYSIS  LOOP 
4HTHE  NUMBER  OF  PEAKS  IS 

3HTYPE  0 - IF  P WAVE  IS  POSITIVE?  TYPE  1 - [ - 3 WAV- 
IVE 

&HWHICH  nL0CK  OF  OATA  FIRST  GRAPHED? 

4 WHO W MANY  BLOCKS  OF  OATA  BETWEEN  EACH  GRA=H? 


Fig.  23.  This  part  of  the  computer  program  sets  parameters 


54 


til 


areas,  form  integer  variables,  rewind  the  data  storage 
device,  and  ask  for  initial  parameters.  In  this  program, 

TAPE  2 specified  in  the  program  identification  statement, 
is  used  as  data  input.  The  output  of  the  program  is  the 
use  of  a plotting  routine  as  specified  in  this  statement. 

The  specification  of  common  storage  was  used  because  of  the 
multiple  use  of  subroutines,  and  the  large  use  of  common 
variables  is  both  the  calling  and  called  routines.  The  P-P 
interval  matrix  (PP}  was  specified  as  an  integer  because 
all  data  values  are  in  integer  form.  The  use  of  REWIND  2 
command  facilitates  using  the  program  many  times  during  a 
single  session. 

The  interactive  part  of  the  program  sets  the  initial 
parameters.  The  variable  IBAUD  is  used  to  set  the  correct 
baud  rate  of  the  CRT  displays  if  used  to  graph  the  output. 
Variable  N is  the  number  of  peaks  in  one  heart  beat  on  the 
ECG  and  is  used  to  decide  where  the  P wave  is  in  relation  to 
the  R wave.  Variable  Nl  will  define  whether  the  P wave  is 
positive  or  negative.  This  is  entered  from  the  ECG  data. 

Nl  is  used  to  decide  which  routine  in  subroutine  READ  is 
used  to  put  the  data  into  an  integer  form.  The  variable  MA 
specifies  the  first  data  block  to  be  graphed.  Variable  MB 
specifies  the  number  of  data  blocks  between  graphing  outputs. 
These  last  two  variables  are  used  to  skip  to  the  data  of 
interest  and  to  control  how  often  the  data  is  graphed.  This 
is  particularly  important  when  computer  time  is  tightly 
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controlled,  as  well  as  allowing  specific  pieces  of  data  output 
t-  be  found  without  all  the  data  being  outputed. 

Read  Data. — Figure  24  shows  the  computer  code  that 
starts  the  data  interpreting  loop  and  reads  the  data  into 
memory.  Loop  counter  K1  specifies  which  data  block  is  being 
read,  analyzed  and  possibly  outputed. 


DO  5 <1=1, 278b 
q*******  READ  DATA 

CALL  REAP 

Fig.  24.  The  code  starting  the  interpreting  loop. 

The  line  CALL  READ  shifts  control  to  subroutine  READ 
(see  Figure  25} . 


*************»**'*i‘'»*-»**»***4',-«  ♦ »**♦■***  If  M»M«<  ► 

SUBROUTINE  REAP  - READ  IN  DATA,  GET  DATA  IMH 
POSITIVE  P WAVE  FORM,  AND  C 0 NV  ERT  TO  INTIGE?  FORMAT 

SUBROUTINE  READ 

COMMON  I DAT  A (51?) , N1 , L< , IP ( 3000 ) ,N 
DIMENSION  DATA (512) 

R EAD ( 2 , i 0 ) (DATA (I) ,1  = 1,512) 

10  FORMAT  (10F8.2) 

I F (EOF ( 2 ) • N E . 0 .) GO  TO  9 
I F ( N1 , ED . 1 ) GO  TO  13 
C*******  FOR  POSITIVE  P WAVE 

DO  1 1=1,512 

IF (DATA ( I ) ,LT.0.)DATA(I)=0. 

1 I DATA  ( I)  =PATA  ( I) 

GO  TO  9 

£♦♦*»***  F OR  NEGATIVE  P WAVE 

13  DO  11  1=1,51? 

I DATA  ( I)  aPATA(I) 

I DATA ( I) =IA8S(I0ATA  (I)  ) 

IF(IOATA(I) .LE.10) IDATA  (I)  = 0 

11  CONTINUE 
N = 2 

9 RETURN 

END 


Fig.  25,  Subroutine  READ  reads  a data  block  of 
512  points  into  memory. 
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Subroutine  READ  used  variable  N1  and  generates  the 
integer  data  matrix  IDATA  and  could  modify  the  value  of  N. 

The  floating  point  data  matrix  DATA  is  generated  and  used 
in  this  subroutine  only.  The  data  is  stored  in  an  eight 
character  field  which  contains  two  decimal  places;  there- 
fore, when  it  is  initially  read  into  memory,  a floating 
point  column  matrix  (1  x 512)  must  be  used.  A test  for  an 
end  of  file  mark  is  made.  If  an  end  of  file  is  encountered, 
the  subroutine  is  exited.  Next,  the  variable  Nl  is  used  to 
decide  which  set  of  routines  will  put  the  data  into  an 
integer  form.  If  the  P wave  is  positive  (Nl  = 0] , the 
values  in  the  DATA  matrix  are  compared,  one  at  a time,  to 
zero.  If  the  data  are  greater  than,  or  equal  to  zero,  the 
value  is  unchanged.  If  the  data  are  negative,  the  value  is 
set  equal  to  zero.  This  can  be  done  because  all  the  infor- 
mation needed  to  find  the  P wave  has  a positive  value.  After 
this  is  done,  the  integer  data  matrix,  IDATA,  is  formed  by 
ignoring  the  decimal  point  in  the  DATA  matrix.  Since  the 
relative  amplitude  of  the  data  s the  only  amplitude  infor- 
mation needed,  this  can  be  done  y. „ „nout  losing  information. 

If  all  of  the  512  datum  points  have  been  handled,  control  is 
returned  to  the  main  program. 

If  the  P wave  is  negative,  the  data  are  handled  dif- 
ferently. First,  the  IDATA  matrix  is  formed  by  ignoring  the 
decimal  point  in  the  DATA  matrix.  Then  the  negative  datum 
points  are  made  positive  by  finding  the  absolute  value  of 
each  datum  point.  Each  of  these  values  are  compared  to  a 
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threshold  (IQ  in  this  case),  that  is  between  the  peak  of  the 
P wave  and  the  peak  of  the  noise  or  Q wave,  whichever  is 
higher.  Every. value  at  or  below  the  threshold  is  set  equal 
to  zero.  After  all  the  data  in  IDATA  have  been  compared, 
the  value  of  N is  set  at  two,  since  the  T wave  usually  has 
a smaller  amplitude  in  this  data  than  the  P wave,  and  thus 
was  nearly  always  below  threshold.  Then  control  was  re- 
turned to  the  main  program. 

Initial  Maximum  Points. — After  the  data  have  been  read, 
a decision  is  made  on  which  data  block  is  being  analyzed  (see 
Figure  26) . If  the  data  block  is  the  first  one,  then  the 


c***,**»  IF  NQT  FIRST  BLOCK  OF  OATA  “I'JD  PEMAINING  ^ w^VES 

IF (Ki»GE»  2) GO  TO  102 

Q**  *********«*****************************, ********************** 

C INITIALLY  FIND  P-P  INTERVAL 

Q***«*******»t«|(**  ************************  »*«**«*»******»***»tt* 

NB±2 

NC-500 

KK=2*N 

c*******  FIND  THE  MAXIMUM  POINTS 

CALL  MAX 


Fig.  26.  Parameters  to  find  the  maximum 

point  on  the  wave  forms  are  set. 


program  does  not  skip  the  next  three  sequences  of  steps. 

The  first  datum  point  in  the  search  window  (.NB)  is  set  equal 
to  two  (2) ; then  the  maximum  datum  point  in  the  search  window 
(NC)  is  set  equal  to  5QG.  Next,  the  number  of  peaks  to  be 
found  (KK)  is  set.  With  these  parameters  established,  sub- 
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routine  MAX  is  called  to  find  tiie  peaks  in  the  data  wave 
forms  (see  Figure  27)  . 

Subroutine  MAX  finds  the  maximum  point  of  a wave  form 
using  a change  in  slope  procedure.  A detailed  explanation 
of  the  theory  used  by  this  subroutine  can  be  found  in 
Chapter  3.  This  subroutine  uses  variables  NB,  NC/  and  KK, 
and  data  matrix  IDATA  to  find  the  wave  form  peaks  and  re- 
turns variable  M and  variable  matrices  IMAX  and  IPP.  In- 
itially, LSIGN  (the  direction  of  the  last  slope)  is  set 
equal  to  one  (1) ; and  M (the  number  of  maximum  points  found) 
is  set  equal  to  zero  (0) . The  size  of  the  data  search  window 
is  set  from  NB  to  NC  and  the  loop  searching  the  data  for  the 
maximum  value  is  started.  The  difference  in  amplitude  be- 
tween the  present  datum  point  and  the  last  datum  point  found 
(NA) . If  NA  is  positive,  the  present  slope  variable  (JSIGN) 
is  set  to  one  (1),  LSIGN  is  set  equal  to  JSIGN,  and  the  loop 
is  continued.  If  NA  is  negative,  JSIGN  is  set  equal  to  zero 
(0) ; and  a comparison  is  made  between  JSIGN  and  LSIGN.  If 
they  are  equal,  the  loop  is  continued.  If  the  difference 
between  LSIGN  and  JSIGN  (LSIGN-JSIGN)  is  less  than  or  equal 
to  zero,  the  loop  is  continued.  If  the  difference  is 
positive,  the  slope  has  turned  negative  and  a maximum  point 
has  been  found.  In  this  case  M is  incremented  by  one  (1) 
and  the  amplitude  of  the  last  datum  point  is  entered  into 
the  maximum  amplitude  matrix  (IMAX) . The  number  of  the 
datum  point  maximum  is  entered  into  the  datura  point  maximum 
matrix  (IPP) . If  this  is  the  first  maximum  point,  LSIGN  is 


£**♦■*•*♦*****'»»****#***■**'***•»«■»■*•**'♦■*■*»**#*'  f »*»*  ♦***^-^*»»  (If  #■  » » * « 4 * 

C SU9R0JTINE'  MAX  - PIND  MAXIMUM  POINT  USING  A CHANG  “ 

C IN  SLOPE  PROCEDURE 

0 *♦♦■*4#  •,*****»*,**.  * *fff«f«fff.f<4«fffffffM>(f*.  fffffffMfffMtfM 

SUBROUTINE  MAX 

COMMON  I DATA  (512) ,N1 ,LK , IP<  3000) ,N 
COMMON  /MAX/ISTEP»IMAX (€) ,1 Pp (6 ) , K< , M B, NC , M 
L SIGN= 1 
M = 0 

DO  13  I=NB,NC 

FIND  THE  OIFFEPENCE  IN  AMPLITUDE  BETWEEN  LAST  AND 
0 PRESENT  DATA  PCINTS  - SLOPE 

NAsIDATA(I)-IOATA<I-l) 

C****-*u  IF  SLOPE  IS  POSITIVE 

IF(NA,GT.O)GO  TO  2 

IF  S10PE  IS  NEGATIVE 
IF(NA»Lr ♦ 0) GO  TO  h 

c****+*k  IF  SL0PE.  Is  7EF0 

GO  TO  13 

F OR  POSITIVE  SLOPE 

2 JSIGN=1 
GO  TO  3 

c*******  FGR  NEGATIVE  SLOPE 

J $IGN=  0 

C*******  IF  SL0PE  OOES  NOT  CHANGE,  DR  IF  SLOPE  GOES  -ROM 

C NEGATIVE  TO  POSITIVE?  CONTINJE  ROUTINE 

I F ( JSIGN.EO.LSIGNJGO  TO  3 
IF( (LSISN-JSIGN) ,LE.O)GO  TO  3 

c**»»***  IF  SL0PE  CHANGE  FROM  POSITIVE  TO  NEGATIVE, SET 

C MAXIMUM  POINT 

M = M + 1 

IMAX(M) alOATA(I-l) 

IPP(M)  =1-1 

c, IF  PjRST  MAXIMUM  FOUND, CONT I NJ E PROGRAM 
IF (M.EO, 1) GO  TO  3 

q*******  IF  ANOTHER  MAXIMUM  FOUND  WITHIN  A REGION,  -[NO 

C MAXIMUM  OF  TWO  VALUES 

I F ( (IPP(M)-IPP(M-l) ) ,GT.3)GO  TO  7 
IFdMAX(M)  .GE.IMAX(M-l)  )5,6 

5 IMAX(M-l)  =IMAX  (M) 

I PP  ( M-  1)  = I PP  (M ) 

6 M=M-1 

C* ******  IF  NUMBER  OF  MAXIMUM  POINTS  EDUAl.  TH  NUMBER  NEEDED, 

C RETURN 

7 IF(M.EQ.K*OGO  TO  15 

S ET  PREVIOUS  SLOPE  INOICATO* 

3 L SIGN=  J3IGN 

13  CONTINJE 

15  RETURN 

END 


Fig.  27.  Subroutine  MAX  finds  the  peaks  of  the  wave  forms. 


set  equal  to  JSIGN  and  the  loop  is  continued.  If  M is 
greater  than  one,  the  difference  between  the  present  and 
last  maximum  datum  point  number  is  found.  If  this  difference 
is  greater  than  three  (3}  (an  empirically  determined  value). 


a comparison  between  M and  KK  is  made.  If  they  are  equal, 


the  subroutine  is  exited.  If  they  are  not  equal,  LSIGN  is 
set  equal  to  JSIGN  and  the  loop  is  continued.  If  the  dif- 
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ference  between  maximum  data  point  numbers  is  less  than  or 
equal  to  three,  the  maximum  amplitude  of  the  two  maximum 
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points  is  found  and  set  into  the  IMAX  matrix  at  position  M-l, 
and  the  number  of  that  maximum  point  is  put  into  the  IPP 
matrix  at  position  M-l.  Next,  M is  reset  to  M-l  and  a com- 
parison between  M and  KK  is  made  as  described  above.  If  NA 
is  equal  to  zero  (0) , the  values  of  LSIGN  and  JSIGN  are  un- 
changed and  the  loop  continues.  When  all  of  the  data  points 
between  NB  and  NC  inclusive  have  been  examined  or  M equals 
KK,  control  is  returned  to  the  calling  routine. 

Finding  Initial  P-P  Interval. — This  part  of  the  computer 
program  identifies  the  P waves  in  the  first  two  heart  beats 
and  finds  the  initial  P-P  interval  (see  Figure  28) . The  P 
wave  counter  (KB)  is  set  to  zero  (0)  first.  Next,  the  DO 
LOOP  that  identifies  the  different  wave  forms  is  initiated. 
This  loop  will  be  repeated  twice,  once  for  each  of  two  heart 
beats.  The  DO  LOOP  counter  is  incremented  by  the  number  of 
wave  peaks  (N) . Therefore,  for  three  peaks  K will  equal  one 
(1)  and  four  (4);  for  two  peaks  K will  equal  one  (1)  and 
three  (3),  but  the  next  increment  is  above  four  so  it  will 
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KB=0 

00  7 <sl,c,N 

C*******  FIND  R AND  THEN  P WAVES 

<BsK3+l 

IF (N»ET»  2) GO  TO  1C8 

IF  t IHAX ( K)  « LE « IMAX  tk*l) ) 109 ♦ 110 

109  IF{ IHAX (< +1 ) ,LE. IHAX <Kf2) ) 111,112 

111  I PC  K9) = I PP ( K +1 ) 

GO  TO  7 

110  IF < IMAX (K)  c IE. IMAX  (K*2)  ) im,  115 

114  GO  TO  111 

115  IP(K9) =IPP(Kf2) 

GO  TO  7 

112  IP(K9) =IPP(K) 

GO  TO  7 

10  8 IFUMAXOO  • LE*  IHAX  { <♦  1)  3 116,117 

116  GO  TO  112 

117  I P(  K8)  =1 PP  ( K *13 
7 CONTINUE 

FIND  p.p  INTERVAL 
PP(l)*ta(2) -IP(l) 


Fig.  28.  The  code  identifying  the  initial  P waves 
and  initial  P-P  interval. 


not  repeat  the  loop.  The  first  step  is  to  increment  KB  by 
one,  then  decide  if  N equals  two.  If  N equals  two,  then  the 
maximum  of  only  two  wave  shapes  is  to  be  found.  Therefore, 
control  is  shifted  further  down  the  program.  The  next  13 
steps  are  just  a series  of  comparison  between  maximum  points 
to  find  the  highest  peak.  This  peak  is  then  labeled  the  R 
wave,  and  the  wave  just  before  it  is  labeled  the  P wave.  If 
the  R wave  is  the  first  wave  form,  the  last  wave  form  in  the 
heart  beat  is  labeled  the  P wave.  The  amplitudes  of  the  P 
and  T waves  are  never  compared  because  these  will  vary  de- 
pending on  EGG  lead  placement.  The  data  point  number  that 
corresponds  to  the  P wave  s entered  into  the  P wave  identity 
matrix  IP.  After  the  P waves  in  the  first  two  heart  beats 


have  been  identified,  the  wave  form  identifying  loop  is 
exited.  The  initial  P-P  interval  is  found  by  subtracting 
IP  (2)  from  IP  (11,  This  is  then  entered  into  the  P-P 
interval  matrix  PP  which  gives  the  P-P  interval  in  the 
number  of  data  points  between  P waves.  This  number  can  be 
multiplied  by  0.002  to  find  the  time  in  seconds  between  P 
waves . 

Parameters  to  Find  Remaining  P Waves . — The  first  five 
equations  set  the  parameters  to  find  the  remaining  P waves 
in  the  first  data  block  (see  Figure  29) . The  next  four 


>******* 


SET  PARAMETERS  TO  FINO  REMAINING  P WAVES  IV 
0 FIRST  SLOCK  OF  DATA 

J = 1 
L K=i 
*NPT=6 
ISTEP=I3(2) 

q *******  SET  PARAMETERS  TO  FIND  P WAVES  IN  REMAINING  3. OCXS 

C OF  DATA 

GO  TO  119 

102  IIPsPP  ((.<)•  (513-ISTEP) 

I ST  EP=  1 
KNPT=LK*3 


Fig.  29.  Parameters  to  find  the  remaining  P waves  are  set. 


equations  modify  almost  all  these  parameters  so  the  P waves 
in  subsequent  data  blocks  can  be  found.  The  variable  J is 
now  the  P wave  counter  for  output  purposes.  The  variable 
LK  serve  as  a counter  on  the  P-P  interval  as  well  as  the 
initial  parameter  on  the  next  DO  LOOP.  Variable  KNPT  serves 
as  the  terminal  parameter  on  the  next  DO  LOOP.  Both  of  these 
will  be  changing  during  execution  of  the  program,  so  they 


63 


must  be  variables.  The  number  of  the  datum  point  being 
examined  is  controlled  hy  ISTEP.  The  first  series  of  com- 
mands in  this  section  equal  the  last  entry  in  the  IP  matrix. 


With  this  set  of  commands  completed,  control  is  skipped 
further  down  the  program. 

If  this  is  not  the  first  data  block  being  examined,  the 
second  set  of  equations  is  used  to  set  the  needed  parameters. 
The  first  equation  in  this  set  is  used  to  find  the  expected 
position  of  the  first  P wave  in  the  new  data  block.  Variable 
IIP  stands  for  the  distance  to  the  expected  P wave.  The 
distance  in  this  case  is  the  number  of  datum  points  left  in 
the  P-P  interval,  PP(),  after  subtracting  the  number  of 
points  in  the  last  data  block  since  the  last  P wave,  513- 
ISTEP . This  number  (IIP)  is  used  in  a later  subroutine  as 
the  center  of  a search  area  to  find  the  P wave  peak.  ISTEP 
is  reset  to  one  (1)  so  the  first  data  point  will  be  con- 
sidered. KNPT  is  always  greater  than  LK  by. five  so  the  next 
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DO  LOOP  will  not  be  exited  before  all  the  P waves  are  iden- 
tified. J is  incremented  to  the  last  identified  PP  interval 


plus  one  (1) . 


Find  Remaining  P Waves . — This  section  of  the  program  is 
a DO  LOOP  that  finds  the  P wave  by  calling  the  subroutine 
GREATER  and  then  calculating  the  P-P  interval  (see  Figure  30) . 
The  parameters  for  the  DO  LOOP  are  LK  and  KNPT  which  are 
described  above.  The  call  to  subroutine  GREATER  will  return 
the  location  of  the  P wave. 
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siLL  mi...* 


Q*»****a  PXN0  p WAVES 

119  00  8 KN=LKtKNPT 

CALL  GREATER 
CALL  GREATER 
IF(N0.GT.512)G0  TO  118 
KL=KN*2 
IP(KL)  =1  STEP 
LK=KN+1 

IF(ISTEP.LT.PP(LK-1))G0  TO  120 
PP(LK) =IP(KL)-IP(LK) 

GO  TO  121 

120  PP(LK) =IP(XL)+LB 

121  IF(ISTE3.GE. 512)60  TO  118 
IIPsPP(LK) 

8 CONTINUE 


Fig.  30.  How  the  main  program  finds 
the  P-P  interval. 


Subroutine  GREATER  (see  Figure  31}  uses  variables  IIP 
and  ISTEP  from  the  main  program,  and  the  matrices  IMAX  and 
IPP  from  subroutine  MAX  to  return  the  P wave  location,  ISTEP, 
to  the  main  program.  Variable  ND  is  the  expected  data  point 
of  the  P wave  peak.  NB  is  the  datum  point  that  is  the  mini- 
mum value  of  the  data  search  window.  Variable  NA  is  the 
initial  value  of  NB  and  will  be  used  later  to  set  the  P wave 
location  if  a peak  is  not  found  within  the  search  window. 

The  maximum  value  of  the  data  search  window  is  NC.  NE  is  the 
initial  value  of  NC.  A test  is  performed  to  see  if  the  maxi- 
mum search  window  value  is  outside  the  data  block.  If  it  is, 
the  subroutine  is  exited.  The  number  of  peaks  to  be  found, 
KK,  is  set  at  six  to  insure  that  the  maximum  peak  in  the 
data  search  window  is  found.  Next,  a call  to  subroutine  MAX 
(see  Figure  27}  is  made.  The  information  returned  from  MAX 
is  the  matrices  IMAX  and  IPP  and  variable  M.  If  one  or  mr  ;e 
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0 SUBROUTINE  GREATER  - FINDS  3 WAVES 

£■***+•*»*»**»»»**  **•♦♦* *^.»». ***»♦*♦»#**♦♦*•**»♦♦♦*•  **♦*.*♦**.4  M< 

SUBROUTINE  GREATER 

COMMON  r DATA  (512)  *N1,LK,IP<  7000  tN/GREU /IIP, NO 
COMMON  /MAX/ISTEP,IMAX(6> ,IPP(&) ,KK,NB, NC,M 
C*******  SET  PARAMETERS 

N 0=  1ST  EP  I p 
NB=ISTE3^I IP-2 
NA  = N3 
NCsNB+5 
N£=NC 


IF(NC.GT ,512)60  TO  201 
<K=6 

£•♦*****  CALL  MAX 


206  CALL  MAX 

IF(M.LE.1)G0  TO  203 

£*»»•*««  if  MORE  THAN  ONE  MAX  POINT, 

DO  204  < A = 1 , M 

IFdMAX(i)  .GE.IMAX(KA)  ) GO  TO  204 
IMAX(l)  =IMAX  (<A) 

IPP(1)=IPP(KA) 

204  CONTINUE 

C*******  OECIOE  IF  P WAVE  IS  PRESENT 

2C3  IF (M.EQ« 1) GO  TO  202 
NB=NB- 2 
NC=NC*2 

IF(  (NC-NB)  -GE.  30)  GO  TO  205 
I F ( NC«  LE « v 1 2 ) GO  TO  207 
NCaNE+1 


FIND  HIGHEST  MAX 


IF(NC.EO.E12)GO  TO  205 
2C 7 GO  TO  206 

C*******  IF  NOT  A MAX  POINT,  KEE3  SHE  D~P  INTERVAL 

205  IBsNAta 
I STE°=I3 
GO  TO  201 

C*******  IF  MAX  POINT,  SET  DATA  INDEX 

202  ISTEP=I:>P(  1) 

201  RETURN 
END 


Fig.  31.  Subroutine  GREATER  returns  the 
location  of  the  P wave. 


peaks  have  been  found,  M will  return  with  a non-oero  number. 
A decision  is  made  if  M is  less  than  or  greater  than  one. 

If  M is  greater  than  one,  the  amplitude  of  the  peaks  are 


compared  in  a DO  LOOP  to  find  the  maximum  peak..  The  number 
of  this  datum  point  is  then  identified  as  the  location  of 
the  P wave.  If  M is  one  or  zero,  a decision  is  made  which 

of  the  two  values  it  is.  If  M is  zero,  then  a peak  was  not 

found  and  the  data  search  window  is  enlarged  by  four  data 
points  (an  empirically  determined  value) . NB  is  decreased 
by  two  data  points  while  NC  is  increased  by  two  data  points. 

If  the  data  search  window  (NC-N3)  is  greater  than  or  equal 

to  30,  the  maximum  size,  then  a P wave  is  assumed  not  to  be 

present.  In  this  case  the  position  of  the  expected  P wave, 

IB,  is  found;  ISTEP  is  set  equal  to  IB;  and  the  subroutine 
is  exited.  If  the  window  is  not  greater  than  30,  a check  to 
see  if  NC  is  within  the  data  block  is  made.  If  NC  is  in  the 
data  block,  the  enlarged  data  search  window  is  examined  for 
a peak  as  described  above.  If  NC  is  outside  the  data  block, 
its  value  is  reset  to  NE  plus  one.  Again  a check  is  made  to 
see  if  the  newest  NC  is  inside  the  data  block.  If  NC  is  equal 
to  the  data  block  maximum  point,  512,  the  P wave  peak  is 
assumed  to  be  at  the  data  block  maximum  as  described  above. 

If  NC  does  not  equal  512,  the  search  is  initiated  again  to 
find  the  P wave  peak  as  described  above.  If  M equals  one, 
ISTEP  is  set  equal  to  the  value  in  IPP  and  control  is  re- 
turned to  subroutine  GREATER. 

When  control  is  returned  to  GREATER,  a decision  is  made 
whether  ND  is  outside  the  data  block.  If  it  is,  then  the 
DO  LOOP  is  exited.  If  ND  is  in  the  data  block,  the  P wave 
location  counter,  KL,  is  calculated.  Next  the  matrix  IP  is 
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set  equal  to  the  value  of  ISTEP  found  in  GREATER.  Next  LK 
is  incremented.  Now  a comparison  between  ISTEP  and  the  last 
value  in  the  PP  matrix  is  made.  If  ISTEP  is  less  than  this 
value,  a new  data  block  has  been  examined  for  a P wave  and  a 
different  procedure  must  be  used  to  find  the  P-P  interval. 
This  procedure  is  to  add  the  location  of  the  present  P wave 


peak  to  the  number  of  data  points  between  the  last  P wave 
and  the  end  of  the  last  data  block  (LB) . If  ISTEP  is  not 
less  than  the  last  value  in  PP,  the  new  value  for  PP  is 
found  by  subtracting  the  present  IP  value  from  the  last  IP 
value.  Independent  of  the  ISTEP  value  in  relation  to  the 
last  PP  value,  a check  is  made  to  see  if  ISTEP  is  outside 
the  data  block.  If  it  is  greater  than  or  equal  to  512,  the 
DO  LOOP  is  exited.  If  ISTEP  is  less  than  512,  IPP  is  set 
equal  to  the  last  PP  value,  and  the  DO  LOOP  is  continued. 

Output. — Output  from  the  main  program  is  two  commands 
(see  Figure  32) . The  first  command  is  a decision  whether 
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Q*¥****4  DECIDE  WHICH  DATA  TO  GRAPH 

118  I E{  MA.  NE  # Kl ) GO  TO  6 
Ml  CALL  MARK 

CALL  MA** 


Fig.  32.  Output  is  basically  a 

call  to  subroutine  MARK 


output  is  desired  by  comparing  MA  and  Kl.  If  MA  does  not 
equal  Kl,  then  control  is  shifted  further  down  the  program. 
If  MA  equals  Kl,  then  output  is  initiated  by  a call  to 
subroutine  MARK. 
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Subroutine  MARK  graphs  the  data,  sets  parameters  for 
marking  P waves,  and  lists  the  values  of  the  P-P  intervals 
in  this  data  block  (see  Figure  33}  . This  subroutine  is  de- 
signed for  use  on  the  Textronix  CRT  terminals  using  the 
Plot-10  package  of  instructions.  Calls  to  subroutines  INITT, 
ERASE,  MOVABS,  DRWABS , ANMODE , and  DCURSR  are  fully  explained 
in  the  instruction  books  for  this  package.  This  subroutine 
can  be  changed  to  graph  the  output  on  any  graphing  terminal 
available.  The  first  command  erases  the  CRT  screen  the  first 
time  there  is  any  output.  Then  IK  is  set  to  47  which  causes 
the  screen  to  be  erased  on  the  next  command  only  in  subse- 
quent outputs.  The  screen  is  prepared  for  data  by  calling 
subroutine  AXIS. 

Subroutine  AXIS  draws  a picture  frame,  a set  of  axes, 
and  the  tic  marks  and  labels  the  axes  and  major  tic  marks 
(see  Figure  34) . The  graphing  subroutines  in  this  subroutine 
are  the  same  as  those  in  MARK.  The  two  data  entries  are  the 
Hollerith  codes  for  labeling  the  Y-axis  (LYLET)  and  major 
tic  marks  on  the  Y-axis  (NYTIC) . The  first  seven  commands 
draw  a 780  by  983  raster  point  frame  on  the  screen.  The 
next  six  commands  draw  the  X and  Y axes.  The  following 
eight  commands  draw  the  minor  tic  marks  on  the  X-axis.  This 
is  followed  by  a DO  LOOP  that  draws  the  major  tic  marks. 

This  organization  of  commands  is  next  repeated  to  draw  the 
minor  tic  marks  first,  and  then  the  major  tic  marks  on  the 
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Y-axis  (see  Figure  35) . After  this,  the  Y-axis  is  labeled 
by  a DO  LOOP  that  prints  one  data  character  in  LYLET  at  a 
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£ ♦♦+*♦+♦***  * *******  * 


3 SUBROUTINE  HARK  - GRAPHS  3A  T A , S FT  PARAMETERS  rDR  MARKING 

C P WAVES » WRITES  P-P  INTERVAL  VALUr 

.++*»**+»*»«»  ♦♦♦♦+**  *****  ♦ *♦*♦-*  «*****«»  . v ***•  ♦****♦♦**  »»♦»»»»*<  ♦ * * ** 

SUBROUTINE  MARK  ■ 

COMMON  I DATA  (512) , Ni, LK  , IP ( 300 0 ) , N/MA RK/I nAUD, J »PP < 2 7i 3 ) 
j** *****  INITIALIZE  SCREEN 

IC(IK.NE.47) CALL  INITT(120) 

IK=47 

CALL  ERASE 

q**.»***  CALL  AXIS 

CALL  AXIS 

C +«*«.+**  PLOT  OATA 

IX*1Q0 
I Y= 150 

CALL  MOVABS (IX  »IY) 

00  300  1=1,512 
NY=IDATA (I) 

CALL  0RWABS<IX+I,IY>NY*2) 

Qmm*  DECIDE  WHETHER  TO  MARK  P WAVES 

IF (IP ( J) .NE.I) GO  TO  300 
LP=200 

LT=IP(J)  >100 

Z*******  CALL  PMARK 

CALL  PMARK (LT» LP) 

J=J  + 1 

300  CONTINUE 
IX  = 15 
IY=600 

CALL  HOVABS (IX  ,IY) 

CALL  MOVA BS ( IX ,IY) 

CALL  ANMODE 

[s***.***#.  WRITE  VALUE  OF  P-P  INTERVAL  I i NUMBER  OF  OATA 


c POINTS.  MULTIPLY  VALUE  BY  .0)2  FOP.  TIME  BETWEEN 

C P WAVES. 

00  301  K=J,LK 
WRITE  90  0 , PP  (K ) 

900  FORMAT  (30X,  16HP-P  INTERVAL  IS  ,13). 

301  CONTINUE 

PAUSE  TO  STUDY  PLOT 
CALL  MOVABS ( 15  ,600) 

CALL  OCURSR(ICH,IX,IY) 

CALL  ERASE 

+ PREPARE  SCREEN  FOR  REMAININ3  3ROORAM 

CALL  MOVABS  (0,780) 

CALL  ANMODE 

303  RETURN 
END 


* *<**  **** 


Fig.  33.  Subroutine  MARK  draws  a graph  of  the  data. 
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SUBROUTINE  AXIS  - DRAW  °ICVJRS  F^AHE,  AXES,  TIC  MA*<$ 
ANO  WRITE  LAPELS 

*«M<4*444»4f>M44«44*4«44«444444M44»*»>  M«4t  >H  ( M4»>  > ) 

SUBROUTINE  AXIS 

£♦**-*,*  0ATA  T0  label  Y-AXIS 

0 ;TA  LYlET/1HA,1HM, 1HP,  lHL,iHItiMT,HU,  IHO.lHE, 

11 H ,7  H (MV/40  ) / 

DATA  NYTIC/2H56,  2H<,2,2H28,2H14/ 

-,***+*,  0RAW  PICTURE  FRAME 

1 X=  15 
I Y=75 

CALL  M OV A 8S ( IX , I Y) 

CALL  0RWABS(IX,IY4-655) 

CALL  DRWARS<IX+1008,IY»655) 

CALL  DRWA9S(IX*1C08,IY} 

CALL  DRW ABS ( IX  ,IY) 

2**>*+  + * DRAW  axis 

I Xs 10  0 
I Ys 150 

CALL  MOVA  PS ( IX ,IY3 
CALL  D RW6 8S ( IX  *8 23 , 1 Y ) 

CALL  MOVARS(IX,IY) 

CALL  DRW4BS<IX,IY+53Q) 

C**4**m  DRAW  X-AXIS  TIC  MARKS  AND  LABEL  TIC  *ARKS 

I Y= 100 

CALL  MOVA«S <IX,IY) 

DO  1 1*110,923,10 
I Y s 15  0 
IX  = I 

CALL  MOVABS(IX,IY) 

CALL  DRW ABS (IX,IY-5) 

1 CONTINUE 

IY=150 

OO  8 K*?10,923,210 
CALL  M0VA8S ( K, IY) 

CALL  ORW A BS ( K,  IY-15 ) 

8 CONTINUF 


(!♦♦♦♦♦**  DRAW  Y-AXIS  ANC  LABEL  TIC  rARXS 

DC  5 I =160 ,680  ,10 
IX=100 
IY  = I 

P ALL  MOVARS ( IX ,IY) 

CALL  0RWABS(IX-5,IY> 

5 CONTINUE 

00  9 1=290,680,170 
CALL  MQVABS (10  0, L) 

CALL  QRWABS (85, L) 

9 CONTINUE 

c»*»**„*  LABEL  Y-AXIS 

CALL  MGVABS (40 ,530) 

CALL  MOV  AqS ( 40 ,5  30 ) 

CALL  ANMOOE 
00  20  1=1,  10 
WRITE  900,LYLET(I) 

900  FORMAT (5X,Ai) 

20  CONTINUE 

WRITE  904, LYLET (11) 

904  FORMAT  (4X».f 
WRITE  905 

905  FORMAT (jX,t*v/») 

c»**44**  cABEL  X-AXIS  MA JOR  TIC  MAR<S 

WRITE  901 

901  FORMAT (4  OX, 3H4  20 ,24X,3H640 ,24X, 44126  0 ) 
WRITE  906 

906  FORMAT  (1 ( / ) ) 

C4444444  LABEL  X-AXIS 

WRITE  902 

902  FORMAT (53X,11HTIME  (MSEC)) 

C4444444  label  y-axis  tic  marks 

IY=685 

CALL  MOVABS (70 ,IY) 

CALL  MOVABS (70 ,IY) 

CALL  ANMOOE 
DO  21  1=1,4 
WRITE  903, NYTIC( I) 

903  FORMAT (10X,A2, 10  (/) ) 

21  CONTINUE 
RETURN 
ENO 


75. 


Subroutine  AXIS  is  continued. 


time  for  the  first  IQ  entries.  Character  11  is  printed  and 
finally  21  vertical  spaces  are  skipped- to  set  up  for  labeling 
the  X-axis  major  tic  marks.  These  tic  marks  are  labeled,  and 
then  some  vertical  space  is  skipped  to  label  the  X-axis. 

Finally  the  Y-axis  major  tic  marks  are  labeled  by  a DO  LOOP 
that  prints  one  data  character  at  a time  from  NYTIC,  then 
skips  10  vertical  spaces  before  printing  the  next  set  of 
characters . 

After  control  is  returned  from  AXIS , subroutine  MARK 
plots  the  data  and  marks  identified  P waves  by  using  a DO  LOOP. 
The  X position  of  the  CRT  light  beam  is  determined  by  the  DO 
LOOP  control  variable.  The  Y position  is  determined  by  multi- 
plying the  individual  values  in  IDATA  by  two  (2)  and  adding 
this  to  the  initial  value  of  IY.  A comparison  is  made  between 
the  next  value  in  the  IP  matrix  and  the  DO  LOOP  control 
variable.  If  these  two  are  not  equal,  the  DO  LOOP  is  con- 
tinued. If  the  two  values  are  equal,  a P wave  has  been  iden- 
tified at  this  position  and  needs  to  be  marked.  The  X position 
(LT)  and  Y position  (LP)  are  established  and  subroutine  PMARK 
is  called. 

Subroutine  PMARK  is  a short  subroutine  that  simply  marks 
identified  P waves  (see  Figure  36) . This  subroutine  takes 
the  parameters  set  in  MARK  and  draws  a vertical  line  10Q 
raster  units  high  to  identify  the  P wave.  Then  control  is 
returned  to  MARI'. 

After  the  P waves  have  been  marked,  subroutine  MARK  in- 
crements the  output  P wave  counter  (J)  and  continues  the 


Fig.  36.  Subroutine  PMARK  marks  P waves  in  the  output. 
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DO  LOOP.  When  the  loop  has  been  completed,  the  CRT  beam  is 
moved  to  a point  above  the  plot  and  the  terminal  mode  is 
changed  to  alphanumeric.  Then  a DO  LOOP  is  started  that 
prints  the  P-P  interval  for  those  intervals  found  in  that 
data  block.  After  this  loop  is  completed,  the  program  is 
halted  by  a call  to  DCURSR.  This  halt  is  set  so  the  operator 
can  study  the  output  and  make  copies  if  needed.  A typed 
zero  (0)  and  RETURN  continues  the  program.  Then  the  light  beam 
is  moved  to  the  upper  left  hand  corner,  and  control  is  re- 
turned to  the  main  program. 

Ending  the  Main  Program. — After  control  is  returned  from 
subroutine  MARK,  MA  is  incremented  by  MB  (see  Figure  37) . 


MA»MA*M3 

6 LBs512-ISTEP 

£****♦*#  PAUSE  TO  STUDY  OUTPUT 

WRITE*, “TYPE  1 WHEN  READY  TO  CONTINUE. " 
WRITE*, “TYPE  2 IF  WANT  TO  STOP." 

READ*, N2 

IFCN2.E1.2JGO  TO  900 
5 CONTINUE 

900  STOP 
ENO 


Fig.  37.  The  main  program  is  ended 
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Next  the  value  of  LB  is  calculated  by  subtracting  ISTEP 
from  512.  Then  a chance  for  the  operator  to  terminate  the 
program  is  offered.  Directions  to  the  operator  are  pre- 
sented. The  value  entered  (N2)  is  read  and  a decision  using 
this  value  is  made.  If  N2  equals  two  (2) , the  program  goes 
to  STOP.  If  N2  does  not  equal  two,  the  DO  LOOP  controlled 
by  Kl  is  continued. 


Computer  Program 

PROGRAM  FI  NO  (TAPE?,  INPUT,  OUTPUT  ,P10D 

**********•***»*******»***  + *******.******»*****,  ****************** 

THIS  IS  MAIN  PROGRAM  TO  FIND  D WAVES  IN  AN  EGG 

*************** ******«******************»,»«* i*-' 

COMMON  I OAT A (512) , Ni ,LK , IP ( 3000 ) , N/GREAT /IIP ,NQ/MAR</ 1 3 A UO, J 
1, P° (2785) 

COMMON  7MAX/ISTEP,IMAX(e> ,1 PP (6 ) , «, N 9, NC , M 
INTEGER  PP 
REWIND  2 

*******  REAQ  IN  THE  IPAUO  VALUE  FOR  INITT 

TAKE  THE  LISTED  BAUD  RATE  AND  DIVIDE  BY  13 
WRITE*, 32HTHE  BAUD  PATE  DIVIDED  BY  10  IS 
READ*, 13 AUD 

c*******  SET  UP  READ  AND  ANALYSIS  LOOP 

WRITE*, 24HTHE  NUMBER  OF  PEAKS  IS 
READ* , N 

WRITE*, 53HTYPE  C - IF  P WAVE  IS  =»0SITIVE:  TYPE  1 - [ r » WAVE 
IIS  NEGATIVE 
READ*, Nl 

WRITE* ,36HWHICH  BLOCK  OF  DATA  FIRST  GRAPHED? 

READ*,  MA 

WRITE* , 4 4HHOW  MANY  BLOCKS  OF  DATA  BETWEEN  EACH  GRA°  H ? 
REA0*,MA 
DO  5 <1*1,2786 
c****,»*  READ  DATA 

CALL  READ 

c*******  if  NOT  FIRST  BLOCK  OF  DATA  ^INO  REMAINING  > WAVES 

IF(K1*GE»  2) GO  TO  102 

Q***4**»  ********  ***********4*,  ***.f********  *****  *****  ************, 

C INITIALLY  FIND  P-P  INTERVAL 

Q****4************** *********»***+*******»*•«*»*****************», 

NB*2 

NC*500 

<K*2*N 

3*****»tt  FIN0  THE  MAXIMUM  POINTS 

CALL  MAX 
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